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Inner-shell Absorption Lines of Fe VI— Fe XVI: A Many-body 

Perturbation Theory Approach 
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ABSTRACT 

We provide improved atomic calculation of wavelengths, oscillator strengths, 
and autoionization rates relevant to the 2^3 inner-shell transitions of Fe VI- 
XVI, the so-called Fe M-shell unresolved transition array (UTA). A second or¬ 
der many-body perturbation theory is employed to obtain accurate transition 
wavelengths, which are systematically larger than previous theoretical results by 
15-45 mA. For a few transitions of Fe XVI and Fe XV where laboratory measure¬ 
ments exist, our new wavelengths are accurate to within a few mA. Using these 
new calculations, the apparent discrepancy in the velocities between the Fe M- 
shell UTA and other highly ionized absorption lines in the outflow of NGC 3783 
disappears. The oscillator strengths in our new calculation agree well with the 
previous theoretical data, while the new autoionization rates are significantly 
larger, especially for lower charge states. We attribute this discrepancy to the 
missing autoionization channels in the previous calculation. The increased au¬ 
toionization rates may slightly affect the column density analysis of the Fe M-shell 
UTA for sources with high column density and very low turbulent broadening. 
The complete set of atomic data is provided as an electronic table. 

Subject headings: atomic data, galaxies: active, techniques: spectroscopic 


1. Introduction 

Since the first astrophysical detection of 2 —3 innershell absorption lines of Fe VII-XII 
in the X-ray spectrum of IRAS 13349-1-2438 obtained with the Reflection Grating Spec¬ 
trometer (RGS) on board XMM-Newton (Sako 2001), this so-called Fe M-shell unresolved 
transition array (UTA) has been identified in many soft X-ray sources. The UTA is mainly 
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comprised of a cluster of lines originating from 2p-3d transitions of M-shell iron ions that 
are located between 15-17 A. These lines, when properly modeled, provide important infor¬ 
mation on the ionization structure, column density, and outflow kinematics of the absorbing 
materials. Behar et al. (2001) calculated a complete set of atomic data for such modeling 
using the Hebrew University Lawrence Livermore Atomic Code (Bar-Shalom et al. 2001, 
HULLAC,), and provided an abbreviated list of transition wavelengths, oscillator strengths, 
and antoionization rates. That complete data set has also been incorporated in some com¬ 
monly used plasma modeling codes. 

In most observations where the Fe M-shell UTA has been identihed, individnal featnres 
from different ionization states are not resolved, mostly due to signihcant turbulent velocity 
broadening (> 700 kms“^) at the sonrce. Low statistics and insnfficient spectral resolntion 
may also hamper the identihcation of individual charge states. Therefore, one has to depend 
on global photoionized plasma models to constrain the physical conditions of the absorbing 
gas. In global hts, where individual UTAs are unresolved, the exact rest-frame position of 
each UTA is particularly important, especially since current photoionization balance calcula¬ 
tions fail to predict the formation of Fe M-shell ions consistently with highly ionized species 
of other elements (Netzer 2003, 2004; Kraemer et al. 2004). 

Recently, Holczer et al. (2005) have reexamined the Chandra High Energy Transmission 
Grating Spectrometer (HETGS) spectrum of the active galactic nucleus (AGN) of NGG 3783, 
by coadding all 900 ks available observations. The low turbulent velocity of the ontflow of 
NGG 3783, the high spectral resolution of the HETGS, and the very long exposure time 
enabled the individual absorption features from different Fe M-shell ions to be resolved. 
The analysis of Holczer et al. (2005) showed that the outflow velocity associated with the 
Fe M-shell UTA appears to be different from those of other highly-ionized species, namely, 
Fe XVH, O VH, and O VIII. In fact, the data prefer a zero velocity for the gas associated 
with Fe M-shell ions, and a —590 kms“^ outflow velocity for the gas associated with more 
highly ionized ions. Therefore, a two component model was proposed where the Fe M-shell 
ions are not part of the AGN outflow, with an important implication of much lower mass 
loss rate. 

As stated in Holczer et al. (2005), the conclusions of that paper depend strongly on 
the calculated HULLAG wavelengths. The velocities of the Fe M-shell UTA can be distin¬ 
guished from those of other ions only if the HULLAG wavelengths are accnrate to better 
than 30 mA, which corresponds to a velocity of ~590 kms“^ at 16 mA. However, it has 
been shown that the theoretical method employed by HULLAG can generally produce un¬ 
certainties of up to 20-50 mA for lines in the 10-20 A band (Brown et al. 2002; Gu 2005). 
Unfortunately, no laboratory measnrements exist for the wavelengths of relevant ions except 
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for a few strong lines of Fe XVI and Fe XV (Brown et al. 2001). For these few lines, the 
HULLAC wavelengths are indeed underestimated by 20-30 mA. However, the strongest and 
cleanest absorption features in the observed spectrum of NGC 3783 are from Fe IX-XI. If 
the HULLAC wavelengths for these ions are also underestimated by similar amounts, the 
implied outflow velocity would be similar to that derived from Fe XVH and O ions. Holczer 
et al. (2005) argued that errors in the calculated wavelengths would be distributed randomly 
in both directions, and therefore, the systematic shift observed in the data is unlikely a result 
of theoretical wavelength errors. Furthermore, it was shown that the HULLAC errors for 
the strongest 2p-3d transitions (those which determine the UTA centroid) of the Fe L-shell 
ions are generally much smaller than 30 mA (Brown et al. 2002; Holczer et al. 2005). The 
present work, however, shows that these indirect indications concerning the HULLAC ac¬ 
curacy for the Fe M-shell UTA wavelengths were misleading and that the assumption that 
these wavelengths could be used for velocity measurements was premature. 

Recently, Cu (2005) showed that a combined configuration interaction and second order 
many-body perturbation theory (MBPT) is capable of predicting wavelengths of 2 —3 
transitions of Fe XVH-XXIV to within a few mA. The 2 —3 lines of the Fe M-shell 
UTA are similar to those of Fe XVH-XXIV, except for the addition of M-shell spectator 
electrons. Following the unexpected result of Holczer et al. (2005), in this paper, we apply an 
improved version of the MBPT method to the 2 —3 transitions of Fe VI-XVI. The resulting 
wavelengths of the UTA lines are found to be systematically larger than HULLAC ones by 
15-45 mA. We also give improved autoionization rates, which are signihcantly larger than 
the previous HULLAC results. The new data are used to model the NCC 3783 spectrum, 
and the outflow velocity derived from Fe M-shell UTA is found to be consistent with that 
associated with Fe XVH and O VH-VHI ions. 

In §2, we describe the detail of MBPT and its implementation. §3 presents the results of 
the calculation, and their comparison with the previous HULLAC data set. A brief summary 
is given in §4. 


2. Theoretical Method 

The present calculation is based on the Rayleigh-Schrodinger perturbation theory for 
a multi-configurational model space. The brief description of the method presented here 
closely follows Lindgren (1974), and readers interested in a more detailed account should 
consult Lindgren (1974). 
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The method tries to solve the Schrodinger equation 


Hdcb'^u — - Efc ' hfc , 


( 1 ) 


with perturbation expansion. The no-pair Dirac-Coulomb-Breit (DCB) Hamiltonian Hdcb 
for an 7V-electron ionic system can be written as (Sucher 1980), 

hdcb = 



where hd{i) is the Dirac Hamiltonian for one free electron, Z is the nuclear charge, r* is the 
radial coordinate of the electron i, Vij is the distance between the electrons i and j, and 
is the frequency-independent Breit interaction 


Bij =- 

^ 2r,- 


cXi ' o;j 


(ai-rij)(aj-rij) 


r~ 


L • ij 

where a-i is a matrix vector constructed from Pauli spin matrices ai 


cti = 


0 (Ti 
(Ji 0 


( 3 ) 

( 4 ) 


Hdcb is split up into a model Hamiltonian and a perturbation V. In the present 
case, a convenient choice is 


Ho = J2^hd{z) + U{r,)] 

i 



( 5 ) 


where U (r) is a model potential including the screening effects of all electrons, whose appro¬ 
priate choice makes V as small as possible. In many applications, U{r) is taken to be the 
Hartree-Fock potential of a suitable configuration. Such a choice simplihes the many-body 
expansion because many terms in the series vanish exactly. However, Hartree-Fock potential 
is non-local and complicates the generation of the zero-th order radial wavefunctions. In the 
present implementation, U{r) is approximated by a local central potential, and is derived 
from a Dirac-Fock-Slater self-consistent held calculation, which minimizes the weighted mean 
energy of all relevant conhgurations. As long as the pertnrbation potential V is kept small 
enongh, the exact choice for U{r) does not affect the hnal resnlts. 

The eigenfnnctions and eigenvalnes of Hq are easily obtained by forming Slater 
determinants from single-electron wavefnnctions once U{r) is determined. A subset of 
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will define a model space M, and the remaining states belong to the orthogonal space N. 
A projection operator P is dehned for M, which produces a state in the model space when 
operated on an eigenfunction of the full Hamiltonian 


*2 = P'S!,, ( 6 ) 

where is generally a linear combination of the subset of that belong to the model 
space M. 

As described in Lindgren (1974), one can dehne a wave operator, H, that transforms 4'° 
back to 

( 7 ) 

The original Schrodinger equation can then be transformed to 

= Et^l. (8) 


This equation defines an effective Hamiltonian in the model space 


Heff = PHqP + pvn, 


( 9 ) 


whose eigenvalues are the true eigenenergies of the full Hamiltonian. The effective 
Hamiltonian is generally non-hermitian, and the eigenfunctions, 4/° are not necessarily or¬ 
thogonal. It follows from the Schrodinger equation and the definition of H that the wave 
operator satisfies 

[n,Ho] = vn-nvn, (lo) 

which is the starting point of the many-body expansion formula. The first order expansion 
of PVQ can be written as 




< > 


■E 

r£N 






j~/r 


(11) 


in terms of matrix elements in the model space M. The first order effective Hamiltonian is 
then 




i, >= H'Scb 


E 


yiryrj 


r&N J 


^0’ 

r 


( 12 ) 


where E['^1,cb =< ^i\HDCB\^j >, and H**” =< >. By solving the generalized 

eigenvalue problem for the hrst order effective Hamiltonian, one obtains the eigenvalues in 
second order. 
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The energy levels obtained in the present work are based on the solution of this gen¬ 
eralized eigenvalue problem. This is slightly different from the method used in Gu (2005), 
where following Vilkas et ah (1999), the total second order energy is expressed as 

Et = Bf' + El (13) 


where is the eigenvalues of PHdcbP in the model space with the corresponding eigen¬ 
vectors and is the second order correction due to correlations in the space N 




blb° 


yiryrj 


r&N i,j&M 


ki kj j^Q 


E^' 


(14) 


This is in fact a good approximation of the generalized eigenvalue solution for If one 

replaces the eigenvectors 6°^ with the true eigenvectors bki of -ffe/p above expression 
becomes the exact eigenvalues of the first order effective Hamiltonian. The uncertainties 
introduced by using the eigenvectors of PHdcbP instead of are generally of higher 
order, making it a valid approximation. Nonetheless, the more accurate method of solving 
the generalized eigenvalue problem for hfg// used in the present work. 

In the present work, we restrict the calculation to the ground state energies of ls^2/®3/'^ 
(1 < O' < 11) configurations, and the energies for states in the ls^2/^3/'^+^ configuration that 
are connected to the ground state by electric dipole transitions. These are the absorption 
lines that appear in AGN where the ions are predominantly in the ground states. The model 
potential are calculated separately for the ls^2/®3/'^ and ls^2/^3/'^’''^ configurations. The 
model space M consists of these configurations, and the N space consists of all single and 
double excitations from them. Because the N space contains not only bound states, but 
also continuum, a finite basis set method is used to evaluate the perturbation expansion to 
convergence. The basis set is derived by imposing the boundary condition Q{ri,)/P{ri,) = b 
for the single electron radial wavefunctions, where Q{r) and P{r) are the small and large 
components of the Dirac spinor, is the boundary radius, which is chosen to be large 
enough to contain the Is, 2/, and 3/ wavefunction amplitudes, and b is an arbitrary constant 
chosen to be K/2rhC, where k is the relativistic angular momentum quantum number of the 
one-electron orbital, and c is the speed of light. With this choice, the boundary condition 
reduces to P'{rb)/P{ri,) = 0 in the non-relativistic limit. 


Several small corrections to the Hamiltonian are also included in the calculations, 
namely, finite nuclear size, nuclear recoil, vacuum polarization, and electron self-energy. 
These are all taken into account with standard procedures of the atomic structure theory. 


In addition to energies, we also calculate the oscillator strengths and radiative rates for 
2/-3/' transitions that are connected to the ground state of each ion. Because the ls^2/'^3/'^^^ 
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states are highly autoionizing, the autoionization rates dominate the total lifetime of these 
levels. We calculate all possible autoionization channels between and 

configurations. In calculating the radiative and autoionization rates, one may either use the 
eigenvectors of PHdcbP, or the generalized eigenvectors bki of 77^^^. We have verified 
that the differences in the two choices are within ~ 20% for the ions considered in the present 
work, and we adopt the results of the latter method, which is generally not very different 
from that of HULLAC. 

For ions from Fe I-V, we have not carried out the MBPT calculation, since MBPT is ex¬ 
pected to be less accurate for near neutral ions, and also because these ions are more complex. 
To give a complete database for the Fe M-shell UTA, we have made simple configuration 
interaction calculations for them, and present the results along with MBPT calculations for 
Fe VI-XVI. 


3. Results 

In Figure 1, we compare the wavelengths and absorption oscillator strengths of the 
present calculation and those of Behar et ah (2001) for Fe VII-XVI. As can be seen, the 
oscillator strengths in the two calculations agree well with each other except for a few isolated 
weak transitions. However, the wavelengths in the present work are systematically larger 
than the HULLAC results by 15-45 mA for some of the strongest lines. The differences 
between the oscillator strength weighted wavelengths in the two calculations are also shown in 
the figure. We have no simple explanation for why the MBPT wavelengths are systematically 
longer than the configuration-interaction method of HULLAC. 

Table 1 gives the complete list of Fe M-shell UTA transitions with the ground states as 
the lower levels. The database is a combination of MBPT results for Fe VI-XVI and config¬ 
uration interaction results for Fe I-V. The tabulated data consist of wavelengths, oscillator 
strengths, radiative transition rates, and autoionization rates. The autoionization rates are 
sums over all relevant autoionization channels. The radiative transition rates are those be¬ 
tween the ground states and the specific upper levels, i.e., possible radiative decays to excited 
states are not included. The omission of these alternative radiative channels should have 
negligible effects on the total natural widths of the upper levels, since autoionization rates 
dominate the radiative rates in most cases, and for a few lines of Fe XVI and XV, where 
radiative rates are comparable to autoionization rates, decays to the ground states are the 
dominant branches. 

In addition to wavelengths, another major difference between the present and the earlier 



calculations of Behar et al. (2001) is in the autoionization rates. The present autoionization 
rates are signihcantly larger than the previous results, especially for lower charge states. 
The discrepancy is caused by additional autoionization channels that are included here. 
Let us take the ls‘^2s‘^2p^3s‘^3p^3d excited levels of Fe IX as an example. In Behar et al. 
(2001), only the autoionization channels that involve the ejection of the 3d electron, i.e. 
toward ls‘^2s‘^2p^{3s3p)^, were considered. The present addition of channels to levels of 
ls‘^2s‘^2p^{3s3p)^3d gives total autoionization rates that are ~3.7 times higher. The increased 
autoionization rates broaden further the natural line width and may have implications for 
the curve of growth if turbulent broadening is particularly small. On the other hand, if the 
source has signihcant turbulent broadening exceeding hundreds of kms“^, as is the usual 
case for AGN outflows, this effect is unlikely to be observed. 

Holczer et al. (2005) recently reanalyzed the HETGS spectra of NGG 3783 in the Fe M- 
shell UTA region. The HULLAG data of Behar et al. (2001) was used in that analysis. The 
major conclusion was that the outflow velocity of the gas associated with Fe M-shell ions are 
signihcantly less than that derived from the Fe XVII, O VII, and O VIII lines. By comparing 
the deepest observed centroid of each ion and the HULLAG-calculated individual-ion centroid 
for Fe IX-XVI, velocities between —374 and 118 kms“^ were derived. The same comparison, 
but now carried out with the present MBPT data is shown in Table 2. Except for a few 
ions where the identihcation of a unique UTA centroid becomes difficult, the velocities range 
between —450 and —780 kms“^. We have htted a uniform outhow velocity (—590 kms“^) 
model to the same NGG 3783 spectrum of Holczer et al. (2005), using the present atomic 
data for Fe VI-XVI. The results are shown in Figure 2. Inspection of Figure 2 shows that 
no systematic discrepancies between the observed and theoretical absorption features for 
the Fe M-shell UTA can be readily identihed, and therefore, the Fe M-shell and other more 
ionized ions appear to belong to the same kinematic system. 

We have been able to identify one single laboratory measurement of Fe XVI lines against 
which we may check the accuracy of the present calculation. Brown et al. (2001) gave the 
wavelengths of the three strongest 2p-3d transitions identihed in the emission spectrum taken 
at the Lawrence Livermore National Laboratory’s electron beam ion trap. One of the three 
lines is blended with a strong Fe XVII line, and the other two lines are unblended, which 
allowed the accurate determination of their wavelengths. The two lines are 2p^3s{J = 1/2)- 
2pi/22p\i2^s3dz/2{,J = 1/2, 3/2) transitions, and have measured wavelengths of 15.208(4) 
and 15.115(6) A. The present work gives wavelengths of 15.2086 and 15.1127 A respectively, 
which agree very well with the measured values. In the contrary, Behar et al. (2001) gave 
wavelengths of 15.1890 and 15.0762 A respectively, which are 20-40 mA too small. Brown 
et al. (2001) also identihed the strongest Fe XV line, 2p^3s^{J = 0)-2pi/22p3y23s^3(i3/2(T = 
1), in their spectrum, but did not give its measured wavelength. Through inspection of 
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Figure 1 in their paper and private communications, we estimate that its wavelength is 
15.353(5) A, which also agrees reasonably well with the present result of 15.3588 A, but 
37 mA larger than the HULLAC result. Unfortunately, the Fe XV and XVI absorption 
lines are relatively weak in the NGC 3783 spectrum, and Fe XVI lines are severely blended 
with Fe XVII lines. Nevertheless, the excellent agreement between the present results and 
the measurement for these few lines lends credibility to the MBPT treatment of the entire 
iso-nuclear sequence. 


4. Conclusions 

We have developed a second order many-body perturbation theory with multi-reference 
model space. We apply the method to the calculation of wavelengths, oscillator strengths, 
and autoionization rates of Fe M-shell UTA arising from 2p-3d transitions originating from 
the ground state of Fe VI-XVI. The wavelengths obtained in the present work are system¬ 
atically larger than the HULLAC calculation of Behar et al. (2001); the present oscillator 
strengths of the strong absorption lines agree well with the HULLAC results; and the au¬ 
toionization rates of Behar et al. (2001) are found to be missing important autoionization 
channels, especially for lower charge states. Using the present data for Fe VI-XVI, we hnd 
no evidence that the outflow velocity of the gas associated with Fe M-shell ions is different 
from that derived from Fe XVH and O VH-VHI ions, as claimed in the recent analysis of 
Holczer et al. (2005) using the HULLAC data. A complete list of Fe M-shell UTA lines 
are given, which include the present MBPT calculation for Fe VI-XVI, and a simple con- 
hguration interaction calculation for Fe I-V. We recommend that this new database should 
be preferred over the earlier HULLAC calculation of Behar et al. (2001) in the analyses of 
future absorption spectroscopy where Fe M-shell UTA is prominent. 

M.F. Cu and S.M. Kahn acknowledges the partial support of NASA grants NAC5-5419 
and NNCC04CL76C. The research at the Technion was supported by The Israel Science 
Foundation (grant no. 28/03). E. Behar thanks the Stanford group for their hospitality 
during a visit in July 2005. We thank Shai Kaspi for making the 900 ks HETCS integrated 
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Table 1. Wavelengths, oscillator strengths, radiative decay rates, and autoionization rates 
of 2 —3 transitions from the ground states of Fe I-XVI. For Fe VI-XVI, the MBPT 
method is used and for Fe I-VI, the standard conhguration interaction method is used. 


Ion 

Upper 

A (A) 

A’" (s-i) 

A“ (s-i) 

A“ (8-1)*^ 

fij 

Configuration (J)^ 


XVI 

25 

15.2652 

1.41E+13 

2.11E+13 

2.63E+13 

9.88E-01 

2pi/238i/23d5/2(./ = |) 


XVI 

26 

15.2086 

2.38E+13 

1.68E+10 

1.39E+13 

8.27E-01 

2pi/23si/23d3/2(./ = o) 


XVI 

27 

15.1127 

8.92E+12 

4.03E+13 

6.33E+13 

6.11E-01 

2pi/23si/23(i3/2(>/ = 2 ) 


XV 

16 

15.3588 

1.78E+13 

3.76E+13 

4.27E+13 

1.89E+00 

2pi/23d3/2(J = 1) 


XV 

14 

15.5965 

7.64E+12 

4.57E+13 

6.93E+13 

8.36E-01 

2p3/23d5/2(J = 1) 


XV 

18 

14.1828 

3.06E+12 

l.lOE+14 

3.24E+14 

2.77E-01 

28i/23p3/2(J = 1) 


XIV 

51 

15.6339 

6.11E+12 

4.51E+13 

3.58E+13 

4.48E-01 

2pi/23pi/23d3/2(J = |) 


XIV 

55 

15.5555 

1.16E+13 

9.90E+13 

8.53E+13 

4.22E-01 

2pi/23pi/23d3/2(J = 2 ) 


XIV 

53 

15.6098 

5.06E+12 

4.88E+13 

9.43E+13 

3.70E-01 

2pi/23p3/23d3/2(J = |) 


XIII 

31 

15.8898 

4.87E+12 

1.76E+14 

2.14E+14 

5.53E-01 

2P3/23pi/23p3/23d5/2(J = 

1) 

XIII 

33 

15.8495 

4.13E+12 

l.OlE+14 

1.30E+14 

4.67E-01 



XIII 

38 

15.7811 

2.31E+12 

2.14E+14 

1.32E+14 

2.59E-01 

2pi/23pi/23p3/23rf3/2(*^ = 

1) 

XII 

95 

16.0253 

4.06E+12 

1.19E+14 

1.26E+14 

2.35E-01 

^Py2^Pl/23py2^<^5/2iJ = 

2^ 

XII 

97 

16.0151 

4.77E+12 

5.39E+14 

1.04E+14 

1.84E-01 

2pi/23p3/23d3/2(J = |) 


XII 

93 

16.0345 

9.30E+12 

2.08E+14 

l.OlE+14 

1.79E-01 

‘^pI/2^pIi2^(I-5/2{J = |) 


XI 

90 

16.1761 

5.21E+12 

3.08E+14 

l.OOE+14 

2.86E-01 

2pi/23p3/23<t3/2(-^ = 3) 


XI 

89 

16.1788 

8.63E+12 

2.98E+14 

9.94E+13 

2.03E-01 

2p3^23Pl/23p|/234/2(-^ = 

1) 

XI 

74 

16.2578 

3.36E+12 

3.07E+14 

1.35E+14 

1.86E-01 

2pi/23pi/23p|/23<t3/2('^ = 

3) 

X 

39 

16.3581 

8.48E+12 

3.75E+14 

1.42E+14 

3.40E-01 

2pi/23p|/23<^3/2(>.^ = |) 


X 

38 

16.3661 

5.59E+12 

4.39E+14 

1.90E+14 

3.37E-01 

2pi/23p|/23<t3/2(^ = |) 


X 

43 

16.2951 

3.47E+12 

7.50E+14 

1.99E+14 

2.07E-01 

^^3/2^^3/2^^5/2('^ — 2 ) 


IX 

3 

16.5392 

1.29E+13 

6.31E+14 

1.72E+14 

1.59E+00 

2pi/23d3/2(J = 1) 


IX 

2 

16.7894 

6.93E+12 

5.85E+14 

9.61E+13 

8.79E-01 

^f*3/2^^5/2('^ “ 


IX 

1 

16.9464 

8.44E-I-10 

5.57E+14 

6.26E-I-13 

1.09E-02 

^^3/2^'^3/2('^ — 1) 


VIII 

27 

16.6784 

1.03E-I-13 

6.37E+14 

2.16E-I-14 

4.30E-01 

2pi/23dg^2(J = |) 


VIII 

16 

16.9289 

4.72E+12 

6.42E+14 

1.97E+14 

3.04E-01 

2pg/23d3/23(l5/2( J = |) 


VIII 

19 

16.8030 

3.38E+12 

6.01E+14 

1.64E+14 

2.14E-01 

2pi/23d^l2Sd^/2{J = |) 


VII 

46 

16.9404 

4.03E+12 

6.26E+14 

2.42E+14 

1.73E-01 

‘2pi/23d^l23ds/2i-^ = 2) 


VII 

26 

17.1276 

2.52E+12 

7.69E+14 

3.74E+14 

1.55E-01 



VII 

25 

17.1289 

3.45E+12 

6.13E+14 

2.90E+14 

1.52E-01 

^?'3/2^‘^3/2^'^5/2(>^ “ 


VI 

24 

17.2656 

3.80E+12 

7.30E+14 

3.88E+14 

2.55E-01 

^P3/2^‘^3/2^'^5/2 V = I) 


VI 

35 

17.2105 

3.31E+12 

5.93E+14 

3.13E+14 

1.47E-01 

~ 2 ) 


VI 

39 

17.1881 

2.72E+12 

7.06E+14 

2.33E+14 

1.21E-01 

^P3/2^'^3/2^'^5/2 V = 2 ) 



^'Total autoionization rates in the HULLAC calculation of Behar et al. (2001). This column is not present 
in the electronic version of this table, and is only included here to illustrated the difference between the 
present and HULLAC results. 

'^The configurations are given in the jj coupling notation, and closed subshells are omitted to save space. 


Note. — Only three strongest absorption lines per ion from Fe VI-XVI are shown here. The entire table 
is included in the electronic version of the paper as a machine readable file. 







Table 2. Best-fit velocities and column densities for ions detected in the 14.9-17.5 A 
region of the HETGS spectrum of NGC3783. 


Ion 

'^Observed 

(A) 

\ a 

^Rest 

(A) 

-^Model*^ 

(A) 

Outflow Velocity ^ 
(kms~^) 

Ion Column Density 
(10^® cm“^) 

Configuration (J) 


o+^ 

15.144 ± 0.003 '' 

15.176 

15.146 

-631 ± 59 

400 ± 60 

4p(J = 1/2, 3/2) 




15.970 ± 0.005 ® 

16.006 

15.977 

-675 ± 94 


3p{J = 1/2, 3/2) 



0+6 

17.351 ± 0.005 

17.395 

17.361 

-759 ± 86 

110 ± 20 

lsbp(J = 1) 




17.161 ± 0.005 

17.199 

17.165 

-663 ± 87 


ls6p(J = 1) 




17.048 ± 0.005 

17.084 

17.050 

-632± 88 


ls7p(J = 1) 



Fe+16 

14.980 ± 0.003 

15.013 

14.985 

-659 ± 60 

3.0 ± 0.5 

2pi/23(i3/2(J = 1) 




15.231 ± 0.002 f 

15.261 

15.232 

-590 ± 39 


2p3/23d5/2(J = 1) 



Fe+16 

15.231 ± 0.002 ‘ 

15.265 

15.232 

-668 ± 39 

1.6 ± 0.4 

2pi/23si/23d5/2(J = 

-) 

2 > 


Fe+14 

15.322 ± 0.006 

15.359 

15.329 

-723 ± 117 

0.4 ± 0.1 

2pi/23d3/2(J = 1) 



Fe+13 

15.569 ± 0.008 8 

15.634 

15.603 

-1247 ± 154 

1.1 ± 0.4 

2pi/23pi/23't3/2('l = 

2> 


Fe+12 

15.844 ± 0.014 g 

15.890 

15.816 

-868 ± 264 

1.2 ± 0.3 

2p|/23pi/23p3/23d5/2(J = 

1) 

Fe+ll 

15.970 ± 0.005 ® 

16.024 

15.977 

-1011 ± 94 

2.3 ± 1.5 

2p|/23pi/23p2^23«'5/2 

(j = 

2 

Fe+lO 

16.154 ± 0.005 

16.179 

16.150 

-464 ± 93 

3.0 ± 1.2 

2pi/23p|/23'^3/2(>l = 

3) 


Fe+9 

16.329 ± 0.005 

16.360 

16.327 

-568 ± 92 

5.5 ± 1.0 

^Pl/2^Pl/2^'^3/2iJ = 

-) 


Fe+8 

16.496 ± 0.004 

16.539 

16.507 

-780 ± 73 

4.0 ± 0.5 

2pi/23d3/2(J = 1) 



Fe+T h 


16.678 



3.0 ± 1.0 

2pii23dli^{J = §) 



Fe+6 h 


17.129 



2.0 ± 0.7 

2pi^23dy^3d^/2{d = 

2) 


Fe+5 h 


17.266 



1.5 ± 0.6 

^P3/2‘^G/2'''^5/2G = 

2’ 


Fe+4 h 


17.328 



1.5 ± 0.6 

^P3/2''G/2'''^5/2G = 

1) 


Fe+3 h 


17.387 



< 0.8 


2^ 



'^For Fe-M ions, these are centroids of the deepest feature in each ionic spectrum. The configuration 
labels given in the last column also correspond to these deepest features. 

^’Centroid in full multi-ion model (Fig. 2). 

‘^Estimated by {Xobserved — X^est)c/Xnest- Errors reflect 90% confidence intervals. 

■^Blend of 0+^ and Fe+^^ at 15.19 A fitted for two distinct features. 

^Unresolved blend of 0+^ and Fe’*'^^ reflected in the high velocity and large error on the Fe+^^ column 
density. 

^Unresolved blend of Fe’*'^® and Fe+^^. 

®These UTAs have multiple minima making centroid determination difficult. (See Fig. 2.) 

''Indirect identification based on best-fit model; Strongly blended with the high-n lines of 0+®. 
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Wavelength (A) 


Fig. 1.— Absorption spectra of 2 —> 3 transitions from the ground state of Fe VII-XVI 
at an arbitrary column density and turbulent velocity of, respectively, 8x10^® cm“^ and 
100 kms“^. Black lines are the present results and red lines are the HULLAC calculation 
of Behar et ah (2001). The value of AA shown in each panel is the difference between the 
oscillator strength weighted wavelengths of the present and HULLAC calculations. 
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Fig. 2.— Chandra HETGS spectra of NGC 3783 in the Fe M-shell UTA region. Black line 
with error bars (plotted every six data points) are the data, red line is the model with a 
nniform outflow velocity of —590 kms“^ using the present atomic data for Fe VI-XVI. 










































































